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Vast developments in quantum technology have enabled the preparation of quantum states with more than a 
dozen entangled qubits. The full characterization of such systems demands distinct constructions depending on 
their specific type and the purpose of their use. Here we present a method that scales linearly with the number of 
qubits, for characterizing stabilizer states. Our approach allows simultaneous extraction of information about the 
fidelity, the entanglement and the nonlocality of the state and thus is of high practical relevance. We demonstrate 
the efficient applicability of our method by performing an experimental characterization of a photonic four-qubit 
cluster state and three- and four-qubit Greenberger-Horne-Zeilinger states. Our scheme can be directly extended 
to larger-scale quantum information tasks. 


I. INTRODUCTION 

Multiqubit states are a basic resource for present and fu¬ 
ture generations of quantum information science experiments. 
In particular, W-qubit stabilizer (or graph) states have well- 
proved utility for one-way ^antum computation and quan¬ 
tum information processing • As the number of particles 
increases, the system and its properties become significantly 
more complex. In order to manipulate and exploit such en¬ 
tangled systems, it is crucial to certify the generated states 
with respect to the ideal stabilizer states. The importance of 
analyzing these quantum resources has led to a variety of the¬ 
oretical works 115141^. Each of them shows certain features of 
the system, e.g., fidelity, purity, and entanglement robustness, 
by using the stabilizer operators or their generators il. Here 
we present a compact approach which allows us to simulta¬ 
neously test the most important properties of the generated 
graph states using a minimal number of measurements. Our 
method utilizes the multiparty Greenberger-Horne-Zeilinger 
(GHZ) theorem ifTTIl for a characterization of the quantum 
state by constructing a Bell-type inequality. In this work we 
briefly introduce nonclassical structures, deflned as the crit¬ 
ical identity products (IDs; discussed in detail in |[l3, [3]) 
and their practical applications for: generalized proofs of the 
W-qubit GHZ theorem, estimation of the fldelity of a state, 
and detection of multi-party entanglement. In the laboratory 
we experimentally generate a four-qubit cluster state and fully 
analyze it through IDs. We proceed in the same way with ex¬ 
perimental three-qubit and four-qubit GHZ states in order to 
illustrate the general utility of IDs. We show how our method 
relates to other methods. 


II. THEORY 

In the Hilbert space of N qubits, nonclassical structures 
related to entanglement, contextuality, and nonlocality were 
recently introduced 111^ . which enable addressing founda¬ 
tional quantum physics topics as well as the characterization 
of states useful for quantum information applications. The 
so-called identity products are the most elementary of these 
structures within the W-qubit Pauli group and form the con¬ 
stituents of the more elaborate nonclassical structures. 


Def 1. IDs are sets of M mutually commuting observables 
(Oi, with i = 1,..., M) whose combined product is ±I (respec¬ 
tively, positive and negative ID). 

Each ID can be represented as a table IDM^, where each 
row is a different W-qubit observable and each column corre¬ 
sponds to a different qubit (see Eig. [T]). The rows are tensor 
products of single-qubit Pauli observables Oq = Yq^Zq} 
and single-qubit identity Iq . When each Oq appears an even 
number of times in all the columns, we call the full set whole 
ID (IDM^); otherwise, we call it partial ID (IDM^). 

Def 2. An ID is maximally entangled if its observables Oi 
cannot be simultaneously tensor factorized into two or more 
separate IDs. It is furthermore critical if no deletion of obser¬ 
vables and/or qubits from the set can result in a smaller ID. 

This sort of entanglement is deflned for a set of mutually 
commuting observables rather than for a particular state 
vector, which we can think of as the Heisenberg-picture 
deflnition of entanglement (see Appendix A). As we will 
see, this deflnition of entanglement is crucial for irreducible 
proofs of the GHZ theorem. 

Each ID is representative of a complete class of equivalent 
IDs under permutations of columns (qubits), and local trans¬ 
formations of qubits’ coordinate systems. Every complete 
class of critical IDs belongs to one or more speciflc classes of 
maximally entangled stabilizer states 10. 

GHZ theorem 

Any class of ID that is whole, negative, and entangled gives a 
straightforward proof of the GHZ theorem for a speciflc class 
of maximally entangled X-qubit states and, consequently, a 
Bell-type inequality violation. Eollowing the X-qubit Mer- 
min inequality 111^ . several different approaches have been 
developed to study the nonlocality of multiqubit states, partic¬ 
ularly graph states iflbl - ITsIl . In all of these works the inequal¬ 
ity is based on stabilizer operators. Remarkably, any whole 
negative entangled ID allows a proof that is irreducible for a 
speciflc class of states and also a generalization of the original 
GHZ theorem. 

Let us consider a joint eigenstate of a whole negative critical 
ID and independent single-qubit measurements {X, X, Z, /} 
on each party. The negativity of the ID guarantees that the 
overall product of the expectation values of the multiqubit 
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Figure 1: (Color online) Critical IDs are minimal sets of 
mutually commuting A^-qubit observables that relate to 
specific multi-qubit states. Each row [e.g., light yellow (light 
gray) circle in (a)] is a different joint observable, where the 
implied tensor product symbols are omitted for compactness, 
while each column [e.g. dark red (dark gray) circle in (a)] 
corresponds to a different qubit. When each single-qubit 
Pauli observable appears an even number of times in each 
column of a negative ID, the set enables us to prove the GHZ 
theorem. The tables represent (a) a whole negative ID related 
to the three-qubit GHZ state, (b) a whole negative ID related 
to the four-qubit linear cluster state, and (c) a partial positive 
ID related to the four-qubit GHZ state. 


observables should be -1 according to quantum mechanics 
(QM). On the other hand, the wholeness of the ID guaran¬ 
tees that the overall product should be +1 in any local hidden- 
variable theory (LHVT), so we obtain the GHZ contradiction 
El. Figures [TJa) and [lib) show two whole negative IDs for 
the three- and four-qubit cases, respectively. Note that this 
type of ID exists only for > 2 and requires measuring at 
most M = N + 1 observables for a critical ID. Starting from 
N = 5, it is possible to find entangled whole negative IDs 
with M < N + 1, giving the most compact demonstration of 
the GHZ theorem; for example, there exist one ID5^ and two 
distinct ID5^ |[l3l . While the original proofs of the GHZ 
theorem depend on the preparation of a particular state, these 
IDs can show the proof using any state within a particular sub¬ 
space. 

ID Bell inequality 

We construct the Bell-type inequality, defining first the cor¬ 
responding Bell’s parameter a for a given negative IDM^ 
as 

M M M 

a = Y,^iOi= Z Oi- Z Oi, (1) 

i i\Xi = l i\Xi=-l 

where Oi are the observables of the ID and Xi are the eigen¬ 
values of a specific (target) eigenstate of the ID. The expecta¬ 
tion value of a according to QM is {a)QM = M. In LHVTs, 
the eigenvalues of each Oq must belong to a noncontextual 
value assignment, and because of wholeness the total number 
of Oq assigned to the eigenvalue -1 must be even. Given this 
constraint, we obtain an upper bound on the expectation value 
of a in LHVTs according to 

{oi)lhvt < M-2, ( 2 ) 

which we call the ID Bell inequality (see Appendix A for 
more details). 


ID entanglement witness 

Any Bell-type inequality can be used to experimentally ver¬ 
ify the correlations within a multiparty state. For the two- 
qubit case the Bell parameter related to the Clauser-Horne- 
Shimony-Holt (CHSH) inequality ll^ is a widely used quan¬ 
tity to characterize sources of two entangled qubits imm. 
In a similar way the A^-qubit ID-Bell inequality can be used 
to certify sources of multi-qubit entangled states. 

We can construct a set of general witness operators for each 
ID This is done by constructing {a) for any partic¬ 

ular class C of states and maximizing over the entire class to 
obtain 7 c = max|^^ec( 7 /;|a|' 0 ). The ID C witness operator is 
then 


W^^=jcl-a, (3) 

which guarantees that (W^^) > 0 for all states in C, while 
clearly (W^^) < 0 onl y fo r states close to the target state 
(assuming jc < M) This includes the so-called 

entanglement witnesses [j^, by letting C be the set of all 
biseparable states, and, more generally, the multipartite 
Schmidt-number witnesses a, by letting C include non- 
biseparable states with different Schmidt numbers than the 
target state. For these specific classes, we can use an existing 
analytic solution l!^ to put an upper bound on 7 c, Fc, 
as shown in Appendix A. However, using this method, we 
obtain a bound that is based solely on the target state, with no 
advantage of considering one ID within the set of stabilizer 
observables over another. In some cases maximizing 7 c 
directly for a particular ID gives a stronger discrimination 
than using Fc. A general analytic method for performing 
this direct maximization is an open question, but numerical 
methods remain feasible for many cases, such as the ones 
presented below. 

ID fidelity estimation 

The measured value of the ID Bell parameter {a)exp enables 
us to put a lower bound on the fidelity of an experimentally 
prepared state l^/;) with respect to the intended eigenstate \no). 
For a general IDM^ (provided that it contains M - 1 inde¬ 
pendent generators from the stabilizer group), we consider the 
case that \ 2 p) is 3 . pure state expressed in the eigenbasis of the 
ID, 


l-ip) =a\Ko)+Zbi\Ki), ( 4 ) 

i=l 

where \ni) are the V - 1 other eigenstates in the basis, and 
kP + TJ=i = 1- Using {a)exp, we obtain a lower bound 
on the amplitude of |/^o} and, consequently, on the fidelity of 
state lip) (see Appendix A for the derivation): 

\a\^>{{a)e,p-M + 4)/4 = FiD. ( 5 ) 

This can be generalized for mixed states by replacing 
the left side of inequality (O with (|ap} = Ej=iCj|ajp, 
which is the weighted average amplitude of \no) among 
the pure states that make up the density matrix plus noise, 
p = co// 2 ^ + with \ipj) being equal to 
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(HHi and 'Zjlo Cj = 1. In practice the bound can be used to 
certify the preparation of a specific quantum state using only 
a maximum of N + 1 measurement settings, without resorting 
to complete quantum state tomography (QST) iS, which 
requires 3^ measurement settings. 

We also want to emphasize that the critical IDs are non- 
classical structures by definition. Critical whole negative IDs 
combine all the above-mentioned quantum properties at once. 
But even noncritical IDs, partial IDs, and/or positive IDs can 
show one or more quantum aspects of the considered eigenba- 
sis. Specifically, any ID that contains N independent genera¬ 
tors, whether it is critical or not, gives us a lower bound on the 
fidelity and can also be used for entanglement discrimination. 

III. EXPERIMENT AND RESULTS 

We apply the ID method to characterize an experimental 
four-qubit cluster state, related to the ID5^ [FigHIb)], where 
the cluster state is a specific class of graph states IJ] . As a fur¬ 
ther demonstration of the functionality of IDs we also analyze 
the three- and four-qubit GHZ states, using the corresponding 
ID43 [Figina)] and ID5^ [FigHtc)], respectively. In order 
to generate these entangled states we use a photonic setup 
(Fig. [2]) in a so-called railway-crossing configuration. Due to 
its compactness and high stability, this arrangement has been 
proven to be very suitable for several experiments ll26l - [^ . 
The scheme is based on a double spontaneous parametric 
down-conversion process (SPDC), bulk optics, and motorized 
tomographic elements to achieve reliable measurements over 
long periods. Additional half-wave plates (HWPs) allow us to 
to switch from the generation of cluster states to GHZ states. 

Four-qubit linear cluster state 

By aligning to produce |0“} entangled pairs in the forward 
direction and in the backward direction (see Fig. [2] and 
Ref. for details), where = (|00) ± |ll))/\/2, we ob¬ 
tain the state: 

ICiin) = (lOOOO) + 10011} + IllOO) - |llll))/2, (6) 

which is equivalent to the linear cluster state up to local uni- 
taries (LU), specifically up to where H - 

{Z+X)l\f2 is the Hadamard gate. The polarizing beam split¬ 
ters (PBSs) and the two interferometers in the setup, which are 
necessary to select the above four terms of the state, reduce the 
four-fold count rate to 0.33Hz. 

Test of GHZ theorem 

Each of the ID5^ measurements is acquired for 4800 s. We 
obtain {a)exp = 3.24 ± 0.05, which shows a violation of the 
ID Bell inequality by 4.8cr and consequently proves the GHZ 
theorem for a four-qubit entangled state [Figl^a)]. More de¬ 
tailed results are reported in Appendix B. The uncertainty, like 
all others reported below, is due to Poissonian counting statis¬ 
tics and constitutes a lower limit for the errors. 

ID entanglement witness 

In order to certify the cluster state through ID entanglement 
witnesses, one constructs 7c for any general pure quantum 



Figure 2: (Color online) A femtosecond-pulsed UV-laser 
beam passes twice through a /3-barium borate (BBO) crystal, 
producing pairs of polarization-entangled photons. The 
photons are emitted in forward and backward directions and 
are recombined on polarizing beam splitters (PBSs). 
Walk-off effects are compensated using HWPs and half-thick 
BBOs. Additional HWPs set the entangled pairs to a selected 
Bell state. By postselecting fourfold coincidence events we 
obtain the desired cluster state or GHZ state. Polarization 
analysis is implemented with motorized tomographic optic 
components. 


state. From the analytic method ll^ we find that to discrim¬ 
inate against all biseparable states (Bi), as well as the four- 
qubit GHZ and W states, T^] 3 i^GHZ,w} = 3 (which also coin¬ 
cides with aLHVT = 3), while to rule out certain other max¬ 
imally entangled four-qubit states T^qc = 4 The mea¬ 
sured value of {a)exp enables us to obtain a negative value for 
0^{Bi GHZ w}) O^iqc)' some cases we can 

find better (more negative) values of ) for some specific 
classes of states by using numerical maximization of {'ip\a\'ip) 
to put an upper bound on 7^. A detailed analysis is re¬ 
ported in Appendix B. In Fig.[3ja) we show a few results of 7c 
obtained via numerical maximization. We consider product- 
states, the GHZ state \GHZ 4 ), the W state IW4), and also 
different types of cluster states, since the linear cluster {Cun) 
is not fully symmetric under the exchange of qubits. In par¬ 
ticular, exchanging the order of the qubits, we evaluate jc for 
the Z cluster |C\) and the shear cluster |Cx). The analytic 
method gives Fjc'^ =3. 

For four qubits there are an infinite number of entanglement 
classes that are inequivalent to one another under stochas¬ 
tic local operations and classical communication (SLOCC) 
M. All of these classes can be given in terms of a rela¬ 
tively small number of continuous entanglement monotones 
m, but a general classification for more qubits is not known. 
A more comprehensive calculation is required to obtain the 
upper bound, 7c for such states. In any event our results for 
(W^^) certify the four-party entanglement and rule out other 
particular maximally entangled four-qubit states. 

Fidelity estimation 
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Figure 3: (Color online) (a) Measured expectation values for 
the ID5^ (on the left) and results of the maximization of 7 c 
for different four-qubit entangled states (on the right). 

|Cx)=(|oooo)+|oioi)+|ioio)-|iiii ))/2 and 
|c^)=(|oooo)+|oiio)+|iooi)-|iiii ))/2 are reached by exchanging 
the order of qubits in the linear cluster state. In the dashed 
box we report the experimental result of the ID-Bell 
parameter, (b) Measured expectation values for the ID4^ (on 
the left) and results of the maximization of 7 c for three-qubit 
entangled states (on the right). In the dashed box we report 
the experimental result of the ID-Bell parameter. 


Using Eq.( [Ml) for the ID5t and {c^)exp, we estimate 
Fjd = 0.56 ± 0.01. Here we want to point out that the 
stabilizer group of the cluster state contains eight different 
ID5^’s that are equivalent by definition, and thus each of 
them allows for a quantum state estimation. All of these 
sets report similar values of Fjd (see Appendix B for the 
complete data). In order to verify the validity of this bound 
we reconstruct the full density matrix through QST with an 
acquisition time of 600s for each measurement setting. The 
extracted quantum state fidelity is Fqst = 0.629 ± 0.007. 

Three-qubit GHZ state 

Measuring one of the cluster state qubits and performing LU 
transformations, we produce the three-qubit GHZ state: 

\GHZ 3 ) = (| 000 > + |111})/V2. (7) 

In the experiment we project the second qubit from Eq.( 6 ) 
onto the diagonal state |-} = (| 0 } - |l })/\/2 and apply a 
Pauli-X operation and a Hadamard gate on the first qubit 
as postprocessing. The state is characterized by the ID4^. 
We analyze it following the same procedure used for the 
cluster state. The GHZ theorem is proven by a violation of 
the ID Bell inequality of 3.la. The ID Bell parameter is 
{<^)exp = 2.6 ± 0.2. We report the 7 c values for the entan¬ 
glement witness in Fig. Ob), with Tc = 2 for biseparable 
states. The obtained {a)exp is not sufficient to rule out 
the three-qubit IIE 3 ); nevertheless, it can still confirm the 
three-party entanglement of the generated state. The fidelity 
values obtained from the ID and QST are Fjd = 0.64 ± 0.05 


and Fqst = 0.672 ± 0.015. Note the relative error for the 
fidelity bound is higher than that for the cluster case, since 
the data are determined from the tomography measurements 
and so are acquired in less time (600s). See Appendix B for 
detailed data. 

Four-qubit GHZ state 

Aligning the two entangled pairs in the setup (Figl2|) to a 
state and a state, with \ 2 p'^) = (| 01 } + | 10 })/a/ 2 , the four¬ 
fold coincidences correspond to the four-qubit GHZ state up 
to two local unitary Pauli-X operations: 

|Gi7Z4> = (|0000> + |llll))/'y2. (8) 

We experimentally implement these LU transformations by 
using HWPs for the third and fourth qubits of the state. The 
state is described by the ID5^ [see Figd^c)], which is critical 
and partial. This implies that the IDs analysis cannot include 
a proof of the GHZ theorem. However, the ID5^ is still max¬ 
imally entangled. It generates the complete stabilizer group 
of the GHZ state, so it can be exploited to test the fidelity of 
the state and as an entanglement witness. We obtain a bound 
of the fidelity of Fm = 0.71 ±0.01 and reconstruct the ex¬ 
act fidelity via QST with the result Fqst = 0.701 ± 0.008. 
The analytic bound for the ID witness, Fc = 3, and {a)exp = 
3.84 ± 0.05, combine to form a negative (W^^) over the class 
of biseparable states. Additionally, the numerical maximiza¬ 
tion technique reports a maximum of 7 c = 3 for several maxi¬ 
mally entangled four-qubit states (see Appendix B for the nu¬ 
merical results), allowing {a)exp to discriminate these from 
the generated state. 


IV. COMPARISON OF DIFFERENT METHODS 

An interesting question is how IDs compare to other ap¬ 
proaches used for state characterization of multiqubit states 
based on incomplete data. 

Concerning the nonlocality proof, we emphasize that the ID 
Bell inequality is composed of a minimal and irreducible set 
of mutually commuting observables for a specific state. This 
is in contrast to previous works [ssl] where the joint ob¬ 
servables are not maximally entangled, implying that nonlo¬ 
cality could still be proven by preparing a state with fewer 
entangled qubits and using fewer parties. While our nonlo¬ 
cality test does not rule out hybrid hidden-variable models of 
entanglement or nonlocality 11 ^ . it does simultaneously dis¬ 
criminate against less entangled states within the Hilbert space 
formalism, as well as some different maximally entangled N- 
qubit states. The Bell inequality for graph states proposed in 
Ref. (Hi involves the complete stabilizer group (SG), which 
is always maximally entangled but is not as compact as an ID, 
scaling exponentially with N rather than linearly. 

Several witnesses were introduced to discriminate specific en¬ 
tangled states BSIll, providing analytic solutions, which 
require minimal experimental effort. Nevertheless, there was 
no generalization for the whole class of stabilizer states, only 
distinct derivations per subclass. For example. Ref. 11] pro¬ 
poses a reduced witness for A^-qubit cluster (GHZ) states 
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which requires (A^ + 1) measurement settings. The ID wit¬ 
ness requires at most N + 1 measurement settings for every 
stabilizer state, and for many specific cases it needs less than 
N settings (e.g. the IDb^ can be measured with four settings 
and the IDbp with only three). Each of these methods is 
minimal in some particular way, and both are robust against 
noise. An additional method for entanglement discrimination, 
discussed in detail in fl^ . is to select subsets of stabilizer ob¬ 
servables that are optimal for discriminating against a partic¬ 
ular state, although a general method for obtaining these sets 
for N qubits is lacking. Unlike critical IDs, these subsets are 
usually not suitable as general entanglement witnesses, be¬ 
cause they do not simultaneously discriminate against other 
particular states or less-entangled states. Reference also 
gives a general method for discriminating between A^-qubit 
stabilizer states using their complete stabilizer groups, but this 
method scales exponentially. The minimal ID witness sets 
can be simultaneously used to discriminate against particular 
states and, in some cases, also to achieve the optimal discrim¬ 
ination against particular states (as with the four-qubit GHZ 
state using the ID4^ in the Appendix B). 

A fidelity estimation with incomplete data is obtained using 
the SG of the state [H, [H, 1^ . This method, based on 2^ 
measurement settings, still scales exponentially, just like the 
QST. Comparing the QST (from ll^ ) and SG analyses for 
our experimental data in FigH] (first two bars), we see that 
the SG fidelity results in a higher value than the QST fidelity 
for states with noise. The QST approach is considered to un¬ 
derestimate the real value of the fidelity llj^ . whereas the SG 
approach, based on the assumption of an a priori known ideal 
state, might jeopardize the actual applicability of the char¬ 
acterized state if the resulting fidelity overestimates the real 
value. Alternatively, a lower bound of the fidelity can be found 
using the generators of the stabilizer group (GoSG) |H,[3^, 
the above-mentioned witnesses (Wit) 11^ . or the IDs. These 
techniques scale linearly and provide thoroughly fair bounds 
for practical applications. Nevertheless the Wit’s derivation is 
not general for stabilizer states like the ID and the GoSG ap¬ 
proaches are. We analytically compare the last two methods 
in Appendix A, showing the IDs give stronger (equally fair) 
bounds on the fidelity within an experimental environment. 
We calculate the fidelity for the experimentally generated sta¬ 
bilizer states using these estimations and summarize the result 
in Fig.m 

We remark that the real value of the IDs approach is to capture 
all the different quantum features of a state at one time. We 
can exploit this generality to calculate the minimum fidelity 
required for an experimental demonstration of multiqubit non¬ 
locality using IDs. Simply setting {a)exp = {ol)lhvt = M-2 
and inverting expression (O, we obtain {\a\^)nonlocal >1/2. 
This verifies the already-proved limit of 50% fidelity, which 
is necessary for violation of any Bell-type inequality based on 
the GHZ theorem In most cases it is also the bound 

for discriminating less than maximally entangled states. 


■ QST BSG HID ■ GoSG HWit 



\Cun) \GHZ^) \GHZ^) 

Figure 4: (Color online) Comparison of fidelities obtained 
with different methods for the four-qubit linear cluster, 
four-qubit GHZ state, and three-qubit GHZ state. The QST 
(red/first bar for every state) and SG (blue/second bar) 
approaches scale exponentially, while the ID (yellow/third 
bar), GoSG (green/fourth bar), and Wit (purple/fifth bar) 
approaches scale linearly with the number of qubits. Within 
the error bars the IDs set lower bounds, in agreement with the 
QST results. The SG fidelities tend to overestimate the QST 
ones. The GoSG and Wit bounds, like the IDs, are consistent 
with the rest of the methods. Note that FqoSG < 0.5 for the 
four-qubit linear cluster, so it is not sufficient to certify that 
the state can violate a Bell-type inequality. The error bars 
derive from Poissonian statistics and thus correspond to a 
lower limit. 


V. CONCLUSION 

We have reported the characterization of an experimental 
four-qubit cluster state and a three-qubit GHZ state with 
the use of critical whole negative IDs. Our efficient method 
requires only A^ + 1 measurements for an A^-qubit state and 
is of high practical value because it provides simultaneously 
a quantum state fidelity bound, an entanglement witness, 
and a nonlocality proof. For these reasons, IDs provide 
convenient laboratory tests of generated entangled resource 
states and certify that they are eligible for quantum science 
applications. Since the ID’s observables belong to a single 
stabilizer group, they can even be implemented within 
stabilizer-based protocols such as quantum error correction 
and measurement-based quantum computing. 

Entangled IDs, even if they are not critical, whole, or negative, 
can still be used to estimate the fidelity of a multiqubit state 
and to construct witness operators, as we have shown with 
the generated four-qubit GHZ state. Additionally, special 
sets of IDs give rise to irreducible proofs of the A^-qubit 
Kochen-Specker theorem ESI, demonstrating the conflict 
between non-contextual hidden variable theories and QM. All 
of these connections emphasize the fundamental relationships 
between entanglement, contextuality, and nonlocality in 
quantum physics. 

Furthermore, in the sense that these nonclassical phenomena 
are exactly the set of resources we wish to exploit, the full 
family of IDs is also the complete set of elemental resources 
for quantum information processing within the A/'-qubit Pauli 
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group. 
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VII. APPENDIX A: THEORY 
A. Derivation of the ID Bell inequality 

In the following we show how to derive the ID Bell inequal¬ 
ity given in Eq.(2) in the main text. 

We rewrite the ID Bell parameter for a given negative IDM^ 
as 

M M M 

a = ^AiOi= Z Oi- Z Ou (9) 

'i i|Ai = l i|Aj=—1 

where Oi are the joint observables of the IDM^ and 
(i = 1,..., M) are the eigenvalues of the ID eigenstate. If a lo¬ 
cal hidden-variable theory (LHVT) is to agree with quantum 
mechanics (QM), then every term Oi in a must be positive 
overall. This means that each Oi with = 1 in Eq.© must 
contain an even number of single-qubit Pauli observables Oq 
assigned the value -1 and each Oi with A^ = -1 must contain 
an odd number of those. Suppose that there are n terms in 
Si|Ai=i Oi that contain two -1 values each, m terms that con¬ 
tain four -1 values each, I terms with six, etc. Likewise, there 
are r terms in Ei|Ai=-i Oi that contain a single -1 value each, 
s terms that contain three -1 values each, t terms with five, 
etc. We also note that because the ID is negative, the value 
7 = r + s + t + ..., which is the overall number of terms in the 
first summation, is always odd. Using these definitions, we 
can write the total number of -1 values appearing in a as 

T] = {2n + 4m + 6/ + ...) + (r + 3s + 5t + ...) 

= (2n + 4m + 6/ + ...) + (2s + 4t + ...)+ 7. (10) 

In the rightmost side of this equation, it is easy to see that 
the numbers in the parentheses are even, and then because 7 
is odd, T] must also be odd. Because the ID is whole, only 
even values of r] are possible in an LHVT, and this causes at 
least one term Oi in a to be negative. Erom this we obtain an 
upper bound, {a)ihvt ^ M - 2, which is finally our ID Bell 
inequality. 


B. Derivation of the ID fidelity bound 

Eor a general IDM^ (provided that it contains M - 1 in¬ 
dependent generators), we consider first the case that |V^} is a 
pure state expressed in the eigenbasis of the ID, 

y 

\lp) = a\Ko) ( 11 ) 

i=l 

where l^c^) are the U - 1 eigenstates in the basis different from 
l^co) and V is the number of all the possible states in the basis. 
Of course |ap + “ 1- Then, the expectation value of 

a 

y 

{(^)exp ~ 1^1 (^o|<^|^o) ■*" ^ \bi\ (12) 

i=l 

We recall that the maximum value of {a)QM is M. Also, 
because the product of all eigenvalues is fixed for the observ¬ 
ables of an ID, any eigenstate \Ki) of the same ID with dif¬ 
ferent values for A^ necessarily causes at least two terms in 
{a) [Eq|9l to be -1, resulting in a maximum of M - 4 for that 
eigenstate. If we allow the presence of noise, Eq.([T^ becomes 

(a)e.p < M\af + ZiM- 4)\bif = 4|a|" + M - 4, (13) 

i=l 

which we can rewrite as 

|aP>((a)e.p-M + 4)/4. (14) 

This is the experimental lower bound on the probability am¬ 
plitude of \ko) within the experimental state \'ip). It corre¬ 
sponds to a lower bound on the fidelity of a particular state 
for M = + 1 and the fidelity that the state lies within a par¬ 
ticular subspace for M < + 1. 

Next, we generalize the above derivation to the case of mixed 
states. Eor a general convex combination of m pure states plus 
noise, 

j m 

p=co^+(15) 

where J] Cj = 1, we can expand each \ij;j) as in Eq.([TT]) , } = 

%l^o) + and follow the same argument to obtain 

m 

(ae.p)< Eci(4|ai|yA^-4). (16) 

j = l 

Given that we have no experimental access to Cj, we must 
allow the constant term to take its maximum value, and then 
we obtain 

{H")>({ae^p)-M + 4)/4, (17) 

where (|ap) = Ecjjajp is the weighted average amplitude 
of |/^o} among the pure states that make up p and the noise 
component (for which the amplitude of |/t:o} is assumed to be 
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ao = 0). Therefore the most general interpretation of our in¬ 
equality is that it places a lower bound on the average ampli¬ 
tude of |/t:o) within a mixed state p and thus that we have ob¬ 
tained a lower bound on the fidelity of the prepared state. This 
also allows for the possibility that our N qubits are entangled 
with additional ancillary qubits that we do not control, since 
measuring them is then analogous to measuring some convex 
mixture of A^-qubit pure states. 


C. Comparing the fidelity bounds obtained using IDs and 
generators 

Let us now compare the fidelity bounds obtained with our 
ID-based method and the generator-based method (GoSG) of 
Ref. JH]. In that work the authors provide a general equa¬ 
tion for any set of N generators which gives the fidelity to be 
bounded below by 

FGoSG = CZ^n-N + 2)l2, (18) 

n 

while our ID-based method gives a lower bound of 

FiD = (Za^-M + A)lA, (19) 

m 

where Ai are observables and = {Ai) are their experimen¬ 
tally obtained expectation values. 

Their method makes use of the N specific generators of a 
graph state, for which all eigenvalues = 1. Every set of 
N independent generators gives an IDM^ (with M = N -\-l) 
by adding one more observable Am to the set, 

= Am Y1 

n 

with Am being equal to the sign of the resulting ID, such that 

0 < < 1 . 

Putting all of this together we can construct a quantitative 
comparison of our two bounds for the same set of N genera¬ 
tors and the Mth observable needed for our method. 

Fid - FgoSG = (^ a™ - M + 4)/4 - (^ a„ - iV + 2)/2 

m n 

= [(aM-l) + (A^-E«")]/4- (21) 

n 

Clearly the difference vanishes when both methods give per¬ 
fect fidelity. However, in the case that the measurements are 
imperfect, -1 < um - 1 < 0 and 0 < - En ^ If we let 

all of the am take the same average value (call it ao < 1), then 
this reduces to 

FiD-FGoSG = {N-l){l-ao)lA>0, (22) 

which shows that our bound is usually better. Of course in 
practice this will depend on the specific values of a^, and 
indeed, in the bizarre case where um = 0 and = 1, we get 
Fid = 0.75 and = 1, and their bound is actually better by 
1/4. So, generally speaking, the best practice will be to take 


the better of these two bounds for a given set of measured 
values am, and their method gives a better bound when 

N - Y^an <1- aM (23) 

n 

or 

(24) 

n 

where = 1 - is the error of each measurement. Interest¬ 
ingly, it is truly arbitrary which of the observables Am in an 
ID is chosen to be Am, which means we can examine all M 
choice, and take the best of the M + 1 different bounds ob¬ 
tained from the measured set a^. Fjd is better for the case 
when the average errors of the all measurements are compara¬ 
ble, but if the error of any one measurement is worse than all 
the others combined, then FqoSG is the superior bound, effec¬ 
tively allowing us to ignore the one particularly bad measure¬ 
ment. The relative quality of the good and bad measurements 
required to satisfy this condition increases linearly with N, 
and thus it becomes increasingly unlikely that we can throw 
away a measurement in this way. Therefore in a realistic ex¬ 
perimental setting, as N increases, Fj d quickly becomes the 
superior bound. 

D. Derivation of the ID entanglement witness 

Here we give the derivation of the analytic solution for the 
upper bound T^ on yc for ID witness observables. We begin 
by rewriting Eq. Cl as 

{G)exp ^ 4|(/co|^)|2 + M-4, (25) 

where \k,o) is the particular eigenstate whose eigenvalues are 
used to define a for this ID. Next, we let C be the class of all 
possible bipartitions {Bi} of the A^-qubit system. Eollowing 
the derivation in we obtain 

max {a)exp < M - 4 + 4[max{z/^}]^ e (Si, (26) 

\P)€Bi m 

where {um} are the Schmidt coefficients of |/^o) with respect 
to the bipartition Bi . We therefore find that T^ = max/ (3i . In 
many cases the individual /3/ have values lower than T^, so 
this method can be used to discriminate more strongly against 
some bipartitions Bi than others. There is also a more general 
analytic solution for jc that rules out some other nonbisepa- 
rable types of states with different Schmidt numbers. 

As in other cases iHlsI], we can also obtain a relation be¬ 
tween these analytic entanglement witnesses and our measure 
of fidelity of the quantum state: 

Fid = {jc-{'Wc^)-M + 4)/4. (27) 

When \^jj) is another stabilizer state, an upper bound can 
also be determined analytically, as shown in M- Eor our 
purposes, this method works by considering which observ¬ 
ables from the ID and the state’s stabilizer act nontrivially 


on the same qubits. For the ID5^ cases presented here, F^ 
gives a bound equal to or better than that of this method; the 
only exception is the case of using ID 5^ (related to the four- 
qubit GHZ state) to discriminate against the shear cluster |Cx} 
(where it gives 7 = 2 while Tc =3, and the numerical re¬ 
sult 7 = 1 is still better). For the ID4^, Tc = 4 is useless 
because that method maximizes over two terms in a sum in¬ 
dependently, ignoring their mutual constraints. In this case, 
the method of lIT^ can still be applied to analytically obtain 
the 7 c = 0 results in Table |II| (see Appendix B), but for all of 
the other cases in that table, it gives 7 = 4, and the numeri¬ 
cal results are still better. This is partly because their general 
method is tailored to discriminating between graph states with 
connected graphs and neglects less entangled states. 

As indicated above, in many cases we can obtain better 
values for 7 c by directly maximizing over {'ip\a\'ip) numeri¬ 
cally. Obviously, no general solution is known for all possible 
classes of states C, but numerical techniques can be used 
to obtain maxima for many particular cases, allowing us to 
discriminate against them, sometimes quite strongly. 

We should point out that these witness techniques implicitly 
assume the Hilbert space formalism of quantum mechanics. 
A more general type of witness can be constructed that rules 
out any hidden-variable theory without pairwise correlations 
between every pair of qubits in the state . Such witnesses 
require one to measure a set of observables that do not all 
mutually commute, so we cannot obtain this result within any 
stabilizer-based protocol. 

E. Noise tolerance of ID entanglement witness 

As has been done in the other cases 10], we can compute the 
general tolerance of our ID witness observables to white noise. 
To compute the tolerance, we solve Tt{W^^p{pn)) < 0 for 
Pat, where p{pn) = + (1 “PAr)|'0}('0| is the standard 

depolarizing noise channel and |V^} is the state we intend to 
witness. For an IDM^ with M = TV + 1, 


More generally, the eigenbasis of an IDM^ is composed of 
projectors |V^)(V^| of rank r = and the noise tolerance 

is given by 


of entangled states. Therefore we do not find a one-to-one 
correspondence between the classification of locally inequiva¬ 
lent entangled graph states (Schrodinger picture) and the clas¬ 
sification of entangled locally-inequivalent IDs (Heisenberg 
picture). This mismatch leads to the existence of maximally 
entangled subspaces (belonging to IDs) that can contain a con¬ 
tinuum of locally inequivalent states (including several locally 
inequivalent graph states). Indeed, the code spaces already 
employed in quantum error correction are of exactly this type, 
although the general utility of maximally entangled spaces is 
more subtle and interesting. 

To get a sense of the structure that emerges here, we can 
look at the four- and five-qubit cases. For four qubits, there 
are three locally inequivalent cluster states (as discussed in 
the main paper); nevertheless, there exists a critical ID4^ with 
an eigenbasis of rank -2 subspaces that contain all three types 
of cluster states. 

For five qubits, there are four classes of maximally entan¬ 
gled stabilizer states up to local unitaries and reordering of 
qubits. These are the five-qubit GHZ state, cluster state, pen¬ 
tagon state, and one other that we will call the cluster-B state. 

The GHZ stabilizers do not contain any ID5^’s, so the five- 
qubit GHZ-type entanglement does not belong to any maxi¬ 
mally entangled subspaces of IDs. The pentagon and cluster 
state share a common negative ID5^, and thus there is a maxi¬ 
mally entangled two-dimensional subspace that contains both 
of these types of states, and all states in this space provide 
proof of the GHZ theorem. There are also critical ID5^’s that 
are common to the cluster and cluster-B states, but none of 
these are whole and negative; thus while they do define maxi¬ 
mally entangled spaces, they do not provide proof of the GHZ 
theorem. There are also numerous spaces that span locally 
inequivalent versions (permutations of qubits) of a given en¬ 
tangled state, just as in the four-qubit cluster case. 

From the above, we can see that the Bell and GHZ states 
look more or less the same in both the Heisenberg and 
Schrodinger pictures, but the same is not true for the other 
types of states. The other types are simply cardinal states 
within complete maximally entangled subspaces that remain 
intact under local unitary evolution. 


VIII. APPENDIX B: ANALYSIS 
A. Four-qubit linear cluster state 


r{M - 7 ) 
r{M - 7) + 7 


(29) 


These tolerances are valid regardless of what method is used 
to obtain 7 . 


F. Entanglement in the Heisenberg picture 

For > 4, there exist maximally entangled IDs with fewer 
than N independent generators that lie at the intersection of 
the stabilizer groups of multiple locally inequivalent classes 


ID entanglement witness 

We present here the method we used to obtain numerical 
bounds for 7 c for the ID5^ in order to discriminate against 
states other than {Cun). 

We break the analysis into pieces based on each LU- 
inequivalent class of an A^-qubit state. This significantly 
reduces the number of parameters needed to explore the gen¬ 
eral state space. For four qubits, a general pure state has 30 
free parameters. If we begin with a particular entangled state, 
then we can explore the entire entanglement class using only 
LU operations, and this reduces the number of free parame¬ 
ters to at most 12 (which is a significant reduction in terms 
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of computational resources needed to calculate the bounds). 
We use the MATLAB OPTIMIZATION TOOLBOX function 
FMINSEARCH.M to perform the multivariate maximization. 
This function finds local maxima based on an initial guess. 
We therefore proceed with a sort of ad hoc “Monte Carlo” 
maximization technique by making a large number of initial 
guesses and taking the best local maximum from among these 
runs. In order to get convergent results from this method, 
we actually compute an upper bound 7 c < max^ Y.T KOi)!- 
This function has far fewer local maxima in C than {a). 
We report in Table [??] the obtained upper bounds of 7 c for 
several quantum states. We considered a fully separable 
state, |V^i)|V^ 2 }|V^ 3 }|V^ 4}5 product states of two-qubit Bell 
states \^ij) {i^j = {1,2,3,4}), partial separable states, GHZ 
states \GHZ), W states \W), and different types of cluster 
states, |Cx} = (lOOOO) + |0101} + |1010} - |llll})/2 and 
\c^) = (lOOOO) + 10110} + |1001) - |llll))/2. 

Although we lack general numerical results for > 4, 
we conjecture that the the negativity of ( which is 

equivalent to the violation of the ID Bell inequality) can 
happen only with the specific stabilizer state that corresponds 
to a (up to LU transformations) - or states that include it as a 
large enough part of a superposition and/or mixed state. 


Z Z Z I 
X X I Z 
Y I X X 
I Y Y Y 

(a) 


State type 

7c 

State type 

7c 

|V’l)IV’2)|V’3>|V’4> 

1 

|V>3>|Gi^^l24> 

2 

k’l>IV’2) $34> 

2 

|V>4>|G/fZl23> 

2 

|V’l)IV’3) $24> 

2 

\^i)\W2Z4) 

2 

V’l) V’4) $23> 

2 

\tp2)\WlZA) 

2 

V’2} V’3) $14> 

2 

\i’z)\Wx24) 

2 

V’2) V’4) $13> 

2 

\tpA)\Wl 2 z) 

2 

k’3)k’4) ^12> 

2 

\GHZ 12 z 4 ) 

0 

|$12)|$34) 

0 

|H/i234> 

2 

j$13)i^’24) 

0 

|Gx) 

4 

i$14)i$23) 

0 

|Gt) 

4 

\^i)\GHZ2za) 

2 

IGlin) 

4 

|V>2)|GifZl34) 

2 




(b) 


Table II: (a) Representation of ID4^ . (b) Numerical upper 
bounds on jc for ID4^. The analytic method always fails 
for this case (i.e., T^ =4). The method proposed in Ref. |[T^ 
gives all cases with 7 c = 0, but fails (i.e. 7 = 4) for all other 
cases. 


Within the cluster stabilizer group there are 196 different 
entangled IDs, belonging to 8 specific isomorphism classes 
with M = 5 or M = 4 and distinct features. From each of 
these we can obtain an ID fidelity and an ID entanglement 
witness. As an example we show in Table lllal one such pos¬ 
itive partial ID4^. The corresponding ID-witness allows us 
to discriminate much more strongly against some entangled 
states with the numerical maximization method than with the 
analytic solution for the same witness [see Tables ?? and ??]. 
Of particular interest are the cases where 7 c = 0 since we can 
discriminate against these states with perfect noise tolerance: 
any {a)exp > 0 is sufficient. 


State type 

7c 

State type 

ic 

|V’l)IV’2)|V’3>|^/’4>* 

2 

\iP3)\GHZi24) 

3 

V’l) V’2} $34> 

3 

|^/>4>|G/fZl23> 

3 

V’l)k’3) ^24)* 

2 

|V’1>|M^234) 

2.6667 

IV’i)IV’4)|$23>* 

2 

|V’2)|M/i34> 

2.6667 

V’2} V’3) $14>* 

2 

|V’3)|VFi24> 

2.3610 

V’2) ^>4) $13>* 

2 

|V’4)|VFi23> 

2.3610 

V’3) V’4) $12> 

2 

\GHZ 12 z 4 ) 

3 

j$12}|$34) 

3 

|VFi234> 

3 

|$13>|$24)* 

1 

|Gx) 

3 

|$14)|$23>* 

1 

Gt_) 

3 

\iJl)\GHZ2Z4) 

\^2)\GHZi34) 

3 

3 

\Cun} 

5 


Table I: Numerical upper bounds on 7 c (max|^)g( 3 (a)) for 
ID5^. All the quantum states, which differ from the target 
state \Ciin) , have the analytic bound F^ = 3, except for 
particular bipartitions (marked with an asterix) where F^ = 2 . 
In some cases the numerical values result are even lower. 


We show the graphs that generate each of the three LU- 
inequivalent four-qubit cluster states in Fig. [5l The graphs 
in Fig. [3b) and (c) are obtained by exchanging the order of 
qubits in the linear cluster state \Ciin) = (| 0000 ) + | 0011 ) + 
| 1100 )-|llll))/ 2 . 



{^\Giin) (b) |Cx) (c) |(7 ) 


Figure 5: Different types of cluster states. 


Quantum state tomography 

We reconstruct the density matrix of the generated cluster 
state through complete quantum state tomography. The real 
part is shown in Fig.jbJ The components of the imaginary part 
are below 0.047 and are hence not presented here. 

The error is estimated running a 100-cycle Monte Carlo simu¬ 
lation with Poissonian noise added to the experimental counts. 
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Z Z 1 1 
Z I X X 
I Z Y Y 

Y X X Y 

Y X Y X 

(a) 

F-min — 0.56 ± 0.01 

Z Z I I 
I Z X X 
Z I Y Y 
X Y X Y 
X Y Y X 

(d) 

Fmin — 0.53 ± 0.01 


Z Z I I 
Z I Y Y 
I Z X X 

Y X X Y 

Y X Y X 

(b) 

Fmin - 0.56 ± 0.01 

I I z z 

Y Y I Z 
X X Z I 

Y X X Y 
X Y X Y 

(e) 

Fmin - 0.52 ± 0.01 


Z Z I I 
Z I X X 
I Z Y Y 
X Y X Y 
X Y Y X 

(c) 

Fmin - 0.53 ± 0.01 

I I z z 

Y Y Z I 
X X I Z 

Y X X Y 
X Y X Y 

(0 

Fmin - 0.51 ± 0.01 


Figure 6: (Color online) Reconstructed density matrix (real 
part) of the four-qubit cluster state (Fqst = 0.629 ± 0.007). 
The imaginary part is not shown since its components are be¬ 
low 0.05. 


I I Z Z 

Y Y I Z 
X X Z I 

Y X X X 
X Y Y X 

(g) 

=0.55 ±0.01 


I I z z 

Y Y Z I 
X X I Z 

Y X Y X 
X Y Y X 

(h) 

F,,,in = 0.54 ± 0.01 


Figure 7: All eight equivalent ID5^ whose joint eigenstate is 

\Clin). 

Stabilizer group 

The stabilizer group operators and their respective expectation 
values are reported in Table HIH 


Observable 

Expectation value 

Observable 

Expectation value 

ZZII 

0.93 ±0.01 

YYIZ 

-0.65 ±0.02 

IIZZ 

0.78 ±0.02 

YYZI 

-0.65 ±0.02 

ZIXX 

0.61 ±0.02 

XYXY 

0.47 ±0.02 

IZXX 

0.59 ±0.02 

XYYX 

0.52 ±0.02 

IZYY 

-0.58 ±0.02 

YXXY 

0.52 ±0.02 

ZIYY 

-0.58 ±0.02 

YXYX 

0.60 ±0.02 

XXZI 

0.66 ±0.02 

ZZZZ 

0.75 ±0.02 

XXIZ 

0.62 ±0.02 

nil 

1±0.03 


Table III: Measured expectation values for all operators in the 
stabilizer group of {Cun). For FqoSG we used the operators 
ZZII, IIZZ, IZXX, XXZL 


Equivalent IDs 

We show in Fig. [7] the eight equivalent ID5^’s within the sta¬ 
bilizer group of \Ciin). We calculate the relative bounds of 
fidelity for each of these IDs, obtaining results in the range 
{0.51±0.01,0.56±0.01}. 


B. Three-qubit GHZ state 

Stabilizer group 

The stabilizer group operators and their respective expectation 
values are reported in Table |IVl Note that these results are 
extrapolated from the quantum state tomography setting of 
the cluster state and after projection of the second qubit of the 
cluster state onto the state |-} = (|0} - |l))/\/2. 


Observable 

Expectation value 

AAA 

0.81 ±0.07 

YXY 

-0.61 ±0.09 

XYY 

-0.59 ±0.09 

YYX 

-0.54 ±0.10 

ZZI 

0.61 ±0.09 

ZIZ 

-0.64 ±0.09 

IZZ 

0.88 ±0.05 

III 

1±0.12 


Table IV: Measured expectation values for all operators in the 
stabilizer group of \GHZs). The first four values are used to 
obtain a Fjd = 0.64 ± 0.04. ZZIJZZ, and XXX are the 
generators used for FgoSG- 
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Quantum state tomography We present in Fig. [ 8 ] the den¬ 
sity matrix of the experimental three-qubit GHZ state, recon¬ 
structed through complete quantum state tomography. 



010 


101 


000 


Figure 8 : (Color online) Reconstructed density matrix (real 
part) of the three-qubit GHZ state (Fqst = 0.672 ± 0.015). 
The imaginary part has components below 0.07 and is not 
shown. 


State type 

7c 

State type 

7c 

|V’l)IV’2)|V’3>|V’4> 

3 

\^3)\GHZi24) 

3 

IV’l) ^>2) $34> 

2 

|^/>4)|GifZl23> 

3 

V’l) V’3) $24> 

2 

\'lpl)\W234) 

2.3333 

V’l) V’4) $23> 

1 

\i’2)\Wl34) 

2.3333 

V’2} V’3) $14> 

1 

\'lp3)\Wl24) 

2.3333 

V’2) V’4) $13> 

2 

\i’4)\Wl23) 

2.3333 

V’s) ^>4) $12> 

2 

\W 1234 ) 

3 

|$12)|$34) 

3 

\C,in} 

3 

j$13)|^’24) 

3 

|Gx) 

3 

i$14)i$23) 

1 

|Gt> 

1 

\iPi)\GHZ234) 

|V>2}|GifZl34) 

3 

3 

\GHZ 1234 ) 

5 


Table V: Numerical upper bounds on jc for ID 5^. For 
biseparable states, the analytic bound is Fc =3, while in 
some cases the numerical result is lower. 

Quantum state tomography We present in Fig. [9] the den¬ 
sity matrix of the experimental four-qubit GHZ state, recon¬ 
structed through complete quantum state tomography. 



1111 

1100 


1100 


1111 


Figure 9: (Color online) Reconstructed density matrix (real 
part) of the four-qubit GHZ state (Fqst = 0.701 ±0.008). The 
imaginary part is not shown since its components are below 
0.03. 


C. Four-qubit GHZ state 


ID entanglement witness 

We report in Table |V| the numerical values of 7 c for the ID 5^ 
calculated via the same maximization procedure used for the 
four-qubit cluster case. The analytic bound is 7 c = 3 for all 
bipartitions. 


Stabilizer group 

The stabilizer group operators and their respective expectation 
values are reported in Table IVll 


Stabilizer 

Expectation value 

Observable 

Expectation value 

ZZII 

0.87 ±0.02 

YYYY 

0.56 ±0.03 

IIZZ 

0.88 ±0.02 

XXYY 

-0.51 ±0.03 

ZIZI 

0.90 ±0.02 

XYXY 

-0.56 ±0.03 

IZIZ 

0.90 ±0.02 

XYYX 

-0.60 ±0.03 

ZIIZ 

0.85 ±0.02 

YXXY 

-0.48 ±0.03 

IZZI 

0.85 ±0.02 

YXYX 

-0.51 ±0.03 

zzzz 

0.85 ±0.02 

YYXX 

-0.53 ±0.03 

xxxx 

0.54 ±0.03 

nil 

1±0.03 


(a) 


(b) 


Table VI: Measured expectation values for the observables in 
the stabilizer group of \GHZ 4 ). The acquisition time for each 
measurement setting was 4800 s. 
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